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We report precision measurements of the Feynman-s (xp) dependence, and first measurements of 
the transverse momentum (pt) dependence, of transverse single spin asymmetries for the production 
of tt" mesons from polarized proton collisions at y^=200 GeV. The xp dependence of the results are 
in fair agreement with perturbative QCD (pQCD) model calculations that identify orbital motion of 
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quarks and gluons within the proton as the origin of the spin effects. Resuhs for the px dependence 
at fixed XF are not consistent with these same pQCD-based calculations. 

PACS numbers: 



The production of particles with high transverse mo- 
mentum from polarized proton collisions at high energies 
is sensitive to the quark (q) and gluon (17) spin structure 
of the proton. Perturbative QCD (pQCD) calculations 
are used to interpret spin observables when they can ex- 
plain measured cross sections. The goal of measuring 
spin observables is to understand how the proton gets its 
spin from its q, g constituents. 

One challenge to theory has been to understand the az- 
imuthal asymmetry of particles produced in collisions of 
transversely polarized protons, known as analyzing power 
(An) or transverse single spin asymmetry. With verti- 
cal polarization, non-zero Apf corresponds to a left-right 
asymmetry of the produced particles. Sizeable An arc 
not expected in coUinear pQCD at leading twist due to 
the chiral properties of the theory Nonetheless, large 
An are observed for inclusive pion production in p-^ + p 
collisions over a broad range of collision energies (v^) 
0, H, fl S Hi and in semi-inclusive deep inelastic scat- 
tering (SIDIS) from transversely polarized proton tar- 
gets [7|. These observations have prompted extensions 
to pQCD that introduce transverse momentum depen- 
dence (TMD) correlated with the spin degree of freedom. 
For example. An could be generated by spin-correlated 
TMD fragmentation if there is transverse q polarization 
in a transversely polarized proton ("Collins effect") Q. 
This mechanism was considered to be suppressed for 
Pf + p ^ TT + X until recently 0, • Spin-correlated 
TMD distribution functions ("Sivers functions") 111. Il2| 
can explain large A n [13| . These functions describe par- 
ton orbital motion within the proton, and so are impor- 
tant for understanding the structure of the proton. 

Although Sivers functions are extracted from SIDIS 
results, there is no proof [lj| that they factorize in pQCD 
calculations ofp|+p— >7r + X. A factorized framework 
involving twist-3 qg correlators has been introduced [l5l | 
and has successfully described [l^ previous An results 
forp|+p TT + X. Of relevance to both approaches 
is a transverse momentum (/cy) that is integrated over in 
inclusive processes. This kx is intrinsic parton motion 
in the Sivers functions and its average is related to the 
inverse proton radius. Large kx is where qg correlators 
are expected to provide a robust framework. Small fc^ 
is where Sivers functions are expected to be applicable. 
Intermediate kx values yield the same results in the two 
approaches, because moments of the Sivers functions are 
found to be related to the qg correlators [13, ■ 

Both theoretical frameworks [H, predict that An 
will increase as the longitudinal momentum (pl) of 
the pion increases, usually given by Feynman-x, xp ~ 
'2pl/\/s- Both frameworks predict that, at fixed xf. An 



will decrease with increasing transverse momentum {pr), 
for pT >1.2 GeV/c. 

Analyzing powers in the hadroproduction of pions have 
been measured before, and typically show a strong in- 
crease as xf increases 

BBSEli- Virtually no previ- 
ous experimental results exist for the dependence of An 
on pt at fixed xf- For ^/s < 20 GeV, the cross sections 
are at least ten times larger than pQC D calculations for 
Xf values where Ajv is sizeable [l9j. This led to the 
suggestion that beam fragmentation, the dissociation of 
the polarized proton by the unpolarized target, was re- 
sponsible for the spin effects, and the expectation that at 
sufficiently large pr these spin effects would vanish. At 
-^5=200 GeV, inclusive tt cross sections at central and 
forward rapidity are found to be in agreement with pQCD 
calculations above pt ^ 2 GeV/c, and are included with 
world data for tt production from e~'"e~ collisions, SIDIS, 
and other p+p collider results in a global analysis of frag- 
mentation functions (20j . Am that increase with xp are 
found at = 200 GeV [1, Im, [111 , but both precision 
measurements and the determination of the dependence 
on Pt have, until now, been missing. 

In this Letter, we report precision measurements of 
the X F dependence and first measurements of the pt de- 
pendence of An at fixed xp for p^ + p tt° + X at 
^/s = 200 GeV. The experiment has been performed at 
the Solenoidal Tracker at RHIC (STAR) ^ at the Rel- 
ativistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Laboratory. The experiment was performed using 
vertically polarized colliding beams. Asymmetries are 
formed from yields measured with left-right symmetri- 
cal detectors, tagged by the polarization direction of one 
beam and summing over the polarization of the other 
beam. Positive Xi? is probed by considering polarization 
of the beam heading towards the detectors and negative 
Xp is probed by considering polarization of the beam 
heading away from the detectors. 

Measurements were carried out with a modular elec- 
tromagnetic calorimeter, known as the Forward Pion De- 
tector (FPD), positioned at large pseudorapidity (?; = 
-ln(tan6'/2)). The (?y)=4.0 results, and some {t])=3.7 
results, reported here were obtained in the 2003 (2005) 
run having integrated luminosity Li„(=0.25 pb~^ (0.1 
pb^^) and average beam polarization ~35% (50%). 
(77) =3.3 and most of the (77) =3.7 measurements were per- 
formed in the 2006 run, which resulted in Li„(=6.8 pb~^ 
with Pb -55%. In the 2006 run, 111 of the 120 possible 
bunches of both RHIC rings, called "Blue" and "Yellow", 
were filled with protons having predetermined patterns 
of polarization signs. The unfilled 9 bunches are sequen- 
tial and correspond to the abort gap needed to eject the 



4 




FIG. 1: Correlation between pion longitudinal momentum 
scaled by ^fsjl (xf) and transverse momentum (pt) for all 
events. Bins in used in Figs. [2] and |4] are indicated by the 
vertical lines. There is a strong correlation between xp and 
Pt at a single pseudorapidity ({??)). 



FIG. 2: Analyzing powers in xf bins (see Fig. [ij at two dif- 
ferent (rj). Statistical errors are indicated for each point. Sys- 
tematic errors are given by the shaded band, excluding nor- 
malization uncertainty. The calculations are described in the 
text. The inset shows examples of the spin-sorted invariant 
mass distributions. The vertical lines mark the tt" mass. 



stored beams. Pt was measured every 3 hours during 
RHIC stores by a polarimeter that detected recoil carbon 
ions produced in elastic scattering of protons from car- 
bon ribbon targets inserted into the beams. The effective 
An of this polarimeter was determined from px + px elas- 
tic scattering from a polarized gas jet target 2J| thereby 
determining Pb = 55.0 ± 2.6% (56.0 ± 2.6%) for the Blue 
(Yellow) beam in the 2006 run 

The FPD comprises four modules, each containing a 
matrix of lead glass (PbGl) cells of dimension 3.8 cm x 
3.8 cm X 18 radiation lengths. Pairs of modules were 
positioned symmetrically left (L) and right (R) of the 
beamline in both directions, at a distance of ~750 cm 
from the interaction point [2l|. The modules facing the 
Yellow (Blue) beam are square matrices of 7 x 7 (6 x 6) 
PbGl cells. Data from all FPD cells were encoded for 
each bunch crossing, but only recorded when the summed 
energy from any module crossed a preset threshold. 



Neutral pions are reconstructed via the decay tt" 



77- 



The offline event analysis included conversion of the data 
to energy for each cell, formation of clusters and recon- 
struction of photons using a fit with the function that 
parameterizes the average transverse profile of electro- 
magnetic showers. Collision events were identified by re- 
quiring a coincidence between the east and west STAR 
beam-beam counters (BBC), as used for cross section 
measurements [23|. Events were selected when two re- 
constructed photons were contained in a fiducial volume, 
whose boundary excludes a region of width 1/2 cell at 
the module edges. Detector calibration was determined 
from the tt" peak position in di-photon invariant mass 
(M^j) distributions. The estimated calibration accuracy 



is 2%. The analysis was validated by checking against full 
PYTHIA/GEANT simulations [26|. The reconstructed 
7r° energy resolution is given by 5Et^/Etj k, 0.16/\/EJJ". 

Due to the limited acceptance there is a strong corre- 
lation between xf and pT for reconstructed tt*^ (Fig. [T]). 
Spin effects in the x f^Pt plane are studied by positioning 
the calorimeters at different transverse distances from the 
beam, maintaining L/R symmetry for pairs of modules. 
Fig.[l]shows loci from (77) =3.3, 3.7 and 4.0. There is over- 
lap between the loci, providing cross checks between the 
measurements. Because the measurements were made at 
a colliding beam facility both xf > and xf < results 
are obtained concurrently. 

Events with 0.08 < M^j < 0.19 GeV/c^ were counted 
separately by spin state from one or the other beam, with 
no condition on the spin state of the second beam, in the 
xp bins shown in Fig. [TJ For each run i, AN,i for each 
bin was then determined by forming a cross ratio 



A 



1 ^Nl^,,Nri^, - ^NLl,^NR^,^ 



(1) 



where A'^/^^^)^^) is the number of events in the 
L (R) module when the beam polarization was up 
(down). Equation ([T]) cancels spin dependent lu- 
minosity differences through second order. Sta- 
tistical errors were approximated by AAjv,i; = 
[Pb^/NL^^^ + NL[,i + Nr^^, + A^i^x,d"^ valid for smaU 
asymmetries. All measurements of Pb for a store were 
averaged and applied to get AN,i for each bin. The run- 
averaged An ± A^jv values are shown in Fig. [2l 

Systematic errors potentially arise from several 
sources. The bunch counter, used for the spin directions. 
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FIG. 3: Analyzing powers versus tt" transverse momentum 
(Pt) for events with scaled tv^' longitudinal momentum I^f] > 
0.4. Errors are as described for Fig. [5] 



identifies events in the abort gaps arising from single- 
beam backgrounds. They account for < 5 x 10^'' of 
the observed yield. Systematic effects from gain varia- 
tions with time are controlled by polarization reversals 
of the stored beam bunches, as demonstrated by ex- 
amples of spin-sorted M^^ for L,R modules in the in- 
set to Fig. [51 Distributions of the significance, Si = 
{A]^ i — Ajv)/A^^r j, are well described by zero mean 
value Gaussian distributions with a equal to unity, as 
expected if the uncertainties are dominated by statistics, 
except near the trigger threshold where larger a is ob- 
served. Systematic errors are estimated from a x A An 
and differences in An associated with tt'^ identification, 
with the largest value chosen. The upper limit on a cor- 
related systematic error, common to all points, arising 
from instrumental effects is 5 An « 4 x 10"**. 

The same pair of modules concurrently measure An 
values consistent with zero for xp < and An that in- 
creases with xp for xp > 0, depending on which beam 
spin is chosen. Null results at xp < are natural since 
a possible gluon Sivers function is probed where the un- 
polarized gluon distribution is large. For xp > 0, a cal- 
culation [isl . [2§| using quark Sivers functions fit (29j to 
SIDIS data [7| best describes our results at (77) =3.3. 
Twist-3 calculations [3l that Ht + p ^ tt + X data 
at ^/s = 20 GeV 01 and preliminary RHIC results from 
the 2003 and 2005 runs at Vs=200 GcV best de- 

scribe the data at (77) =3.7. Both calculations are in fair 
agreement with the variation of with xp. Neither 
calculation describes data at both {rj) . 

Events from modules at different (77) that overlap in the 
xp-pT plane (Fig. [TJ provide consistent results. Hence, 
it is possible to further bin the results not only by 2;^? but 
also by pt- For this analysis, pt is determined from the 
measured energy, the fitted position of the tt*^ within an 
FPD module, and the measured position of the module 
relative to the beam pipe and to the collision vertex. The 
z component of the event vertex uses a coarse time differ- 
ence between the east and west BBC, and is determined 
to ~20 cm resulting in Apt/pt=0.04, where Apr is the 
uncertainty in pT- One method of determining the pt 
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FIG. 4: Analyzing powers versus tt" transverse momentum 
(pt) in fixed xp bins (see Fig. [!}. Errors are as described for 
Fig. [5] The calculations are described in the text. 



dependence (Fig. [3|) was to select events with \xp\ > 0.4. 
An is consistent with zero for xp < —0.4. For xp > 0.4, 
there is a hint of an initial decrease of An with pT, al- 
though the statistical errors are large, since (7y)=4.0 data 
were only obtained in the 2003 and 2005 runs with lim- 
ited integrated luminosity and polarization. For Pt > 1.7 
GcV/c, An tends to increase with pp for xp > 0.4. This 
is contrary to the theoretical expectation that An de- 
creases with Pt- 

The results in Fig. [3] may still reflect small correlations 
between xp and pp for each point, rather than the de- 
pendence of An on pp at fixed xp. To eliminate this cor- 
relation, event selection from Fig. [1] was made in bins of 
Xp, followed by bins in pp. The resulting variation of Am 
with pp is shown in Fig. 21 compared to calculations [13[ 
using a Sivers function fit to pi +p n + X data Q and 
twist-3 calculations For each point, the variation of 
(xp) is smaller than 0.01. There is a clear tendency for 
An to increase with pp, and no significant evidence over 
the measured range for An to decrease with increasing 
Pp, as expected by the calculations. This discrepancy 
may arise from unexpected TMD fragmentation contri- 
butions, xp,pp dependence of the requisite color-charge 
interactions, evolution of the Sivers functions, or from 
process dependence not accounted for by the theory. 

In summary, we have measured the xp and pp de- 
pendence of the analyzing power for forward tt'^ produc- 
tion vn p^ + p collisions at v^=200 GeV in kinematics 
(0.3 <xp < 0.6 and 1.2 < pp < 4.0 GeV/c) that strad- 
dle the region where cross sections are found in agree- 
ment with pQCD calculations. The xp dependence of 
the TT*^ An is in fair agreement with both a collinear 
twist-3 calculation and a calculation assuming factoriza- 
tion that attributes the spin effects to spin-correlated in- 
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trinsic transverse momentum of the quarks within the 
proton. Recent theoretical work interrelates these de- 
scriptions. Both calculations expect the spin effects to 
monotonically decrease with increasing px for px >1.2 
GeV/c. Measurements of the px dependence at fixed xp 
of A]\j are not consistent with these expectations. This 
may reflect the presence of additional mechanisms for 
these spin effects. Future measurements capable of dis- 
entangling TMD fragmentation and distribution function 
contributions to tt^ spin effects, and measurements of Ajsj 
for real and virtual photon production sensitive to only 
Sivers contributions, arc required to definitively estab- 
lish if partonic orbital motion is the correct explanation 
of these effects. 
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